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Abstract 

We introduce 2-3 finger trees, a functional representation of persistent sequences support- 
ing access to the ends in amortized constant time, and concatenation and splitting in time 
logarithmic in the size of the smaller piece. Representations achieving these bounds have 
appeared previously, but 2-3 finger trees are much simpler, as are the operations on them. 
Further, by defining the split operation in a general form, we obtain a general purpose 
data structure that can serve as a sequence, priority queue, search tree, priority search 
queue and more. 



1 Introduction 

Lists are the functional programmer's favourite data type. They are well supported 
by most if not all functional programming languages. However, one could argue 
that this support sometimes leads programmers to use lists when a more general 
sequence type would be more appropriate (Okasaki, 2000). The operations one 
might expect from a sequence abstraction include adding and removing elements at 
both ends (the deque operations), concatenation, insertion and deletion at arbitrary 
points, finding an element satisfying some criterion, and splitting the sequence into 
subsequences based on some property. Many efficient functional implementations of 
subsets of these operations are known, but supporting more operations efficiently 
is difficult. The best known general implementations are very complex, and little 
used. 

This paper introduces functional 2-3 finger trees, a general implementation that 
performs well, but is much simpler than previous implementations with similar 
bounds. The data structure and its many variations are simple enough that we 
are able to give a concise yet complete executable description using the functional 
programming language Haskell (Peyton Jones, 2003). The paper should be acces- 
sible to anyone with a basic knowledge of Haskell and its widely used extension to 
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multiple-parameter type classes (Peyton Jones et at, 1997). Although the structure 
makes essential use of laziness, it is also suitable for strict languages that provide 
a lazy evaluation primitive. 

The rest of the paper is organized as follows. The next section contains some 
preliminary material used throughout the paper. In Section 3 we describe the finger 
tree type, an elegant example of a first-order nested type, and show how it supports 
deque operations and concatenation. In Section 4 we define a variation of the finger 
tree type that supports a splitting operation, abstracting over the annotation used 
to guide splitting. In this way we obtain a general purpose data structure that can 
serve as a sequence, priority queue, search tree, priority search queue and more. 
Moreover, these trees perform well in each role, though specialized implementations 
are usually faster. Section 5 concludes, including a discussion of related work. 

We have included several exercises to support the digestion of the material and 
to discuss further variations and extensions. 

2 Preliminaries 

This section briefly reviews two fundamental concepts, monoids and reductions, 
that are at the heart of the later development. 

2.1 Monoids 

A type with an associative operation and an identity forms a monoid. 

class Monoid a where 
0 :: a 

(©) :: a — > a — > a 

We shall introduce a variety of monoids later on. 

A single type is sometimes the carrier of several different monoids. In this case, 
we will introduce an isomorphic type via a newtype declaration for each monoidal 
structure of interest. 

2.2 Right and left reductions 

Monoids are often used with reductions. A reduction is a function that collapses a 
structure of type / a into a single value of type a (Meertens, 1996): empty subtrees 
are replaced by a constant, such as 0; intermediate results are combined using a 
binary operation, such as A reduction with a monoid yields the same result 
however we nest the applications of '©', but for a reduction with an arbitrary 
constant and binary operation we must specify a particular nesting. If we arrange 
that applications of the combining operation are only nested to the right, or to the 
left, then we obtain skewed reductions, which we single out as a type class. 

class Reduce f where 

reducer :: (a — > b — > b) — > (/ a —> b — > b) 
reducel :: (6 — > a — > b) — > (6 — > / a — > b) 
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Skewed reductions are far more general than reductions: they not only take an 
arbitrary binary operation and an arbitrary element, the arguments of the binary 
operation may even have different types. 

The order of reducer's parameters is arranged as to suggest an alternative reading: 
reducer (<) lifts the operation '-<' to /-structures. As an example, if '-<' adds an 
element to a container, then reducer (<) adds a complete /-structure. 

Skewed reductions need not be written by hand; they can be defined generically 
for arbitrary /-structures (Hinze & Jeuring, 2003). Nonetheless, we shall provide 
all the necessary instances in order to keep the paper self-contained. 

On lists, reductions specialize to folds. 

instance Reduce [] where 

reducer (<) x z = foldr (<) z x — NB. x and z are swapped 
reducel (>) x z = foldl (>-) x z 

Reductions are particularly useful for converting between structures. 

toList :: (Reduce /) => / a — > [a] 
toList s — s :' [] where (:') = reducer (:) 

Think of the structure s as being cons'ed to the empty list. 



3 Simple sequences 

We first develop the definition of 2-3 finger trees, and then show how they efficiently 
implement sequences. As a starting point, consider ordinary 2-3 trees, like the one 
below: 




The structural constraint that all leaves occur at the same level may be expressed 
by defining a non-regular or nested type (Bird & Meertens, 1998), as follows: 

data Tree a = Zero a \ Succ (Tree (Node a)) 
data Node a — Node2 a a \ NodeS a a a 

Values of type Tree a have the form Succ 11 (Zero t) for some n, where t has type 
Node n o, the type of well-formed 2-3 trees of depth n. Note that the internal nodes 
of these trees contain no keys; in this section we shall store all data in the leaves of 
our trees. 

Operations on such trees typically take time logarithmic in the size of the tree, 
but for a sequence implementation we would like to add and remove elements from 
either end in constant time. 
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3.1 Finger trees 

A structure providing efficient access to nodes of a tree near a distinguished location 
is called a finger (Guibas et al., 1977). In an imperative setting, this would be done 
by reversing pointers. In a functional setting, we can transform the structure in a 
manner reminiscent of Huet's "zipper" structure (1997). 

To provide efficient access to the ends of a sequence, we wish to place fingers at 
the left and right ends of this tree. Imagine taking hold of the end nodes of the 
example tree above and lifting them up together. We would obtain a tree that looks 
like this: 




Because all leaves of the original 2-3 tree were at the same depth, the left and right 
spines have the same length, and we can pair corresponding nodes on these spines 
to make a single central spine. Hanging off the sides of these nodes are 2-3 trees, 
whose depth increases as we descend the central spine. The first level contains two 
or three elements on each side, while the others have one or two subtrees. At the 
bottom of the spine, we have either a single 2-3 tree or none, depending on whether 
the old root had a degree of 3 or 2. We can describe this structure as follows: 

data FingerTree a = Empty 
| Single a 

| Deep (Digit a) (FingerTree (Node a)) (Digit a) 

where a digit is a buffer of elements stored left to right (in the picture digits are 
depicted by filled circles), here represented as a list to simplify the presentation: 

type Digit a = [a] 

As noted above, in transformed 2-3 trees these lists have length 2 or 3 at the top 
level and 1 or 2 at lower levels. We shall relax this, allowing between one and four 
subtrees at each level. As we shall see in the next subsection, this relaxation provides 
just enough slack to buffer deque operations efficiently. By the way, the buffer type is 
called Digit because finger trees are a so-called numerical representation (Okasaki, 
1998), a data structure that is modelled after a number system. 

The non-regular definition of the FingerTree type determines the unusual shape 
of these trees, which is the key to their performance. The top level contains elements 
of type a, the next of type Node a, and so on: the nth level contains elements of type 
Node n a, namely 2-3 trees of depth n. Thus a sequence of n elements is represented 
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by a tree of depth 8(logn). Moreover, an element d positions from the nearest end 
is stored at a depth of 0(logd) in the tree. 

Before moving on, we instantiate the generic definition of reduction to these 
types. Reduction of nodes is trivial: 

instance Reduce Node where 

reducer (<) (Node2 a b) z = a < (b < z) 
reducer (<) (NodeS a b c) z = a < (b -< (c < z)) 
reducel (>-) z (Node2 b a) = (z >- 6) >- a 
reducel (>-) z (NodeS c b a) = ((z >- c) >- b) >- a 

Reduction of finger trees uses single and double liftings of the binary operation: 

instance Reduce FingerTree where 
reducer (<) Empty z = z 

reducer (<) (Single x) z = x < z 

reducer (<) (Deep pr m sf) z = pr -<' (m -c" (sf <' z)) 
where (-<') = reducer (<) 

(<") = reducer (reducer (-<)) 

reducel (>-) z Empty = z 

reducel (>) z (Single x) = z y x 
reducel (>-) z (Deep pr m sf) = ((z >-' pr) >-" m) >-' sf 
where (>-') = reducel (>-) 

(>-") = reducel (reducel (>-)) 



We shall now proceed to show that finger trees make efficient deques, with all op- 
erations taking 0(1) amortized time. This holds even if the structure is used in 
a persistent manner, where several versions of a tree may coexist (Driscoll et al., 
1989). The analysis is essentially the same as that for Okasaki's implicit deques 
(1998). Indeed, the FingerTree type differs from implicit deques primarily in stor- 
ing Nodes instead of pairs, with a corresponding increase in the size of digits. This 
change provides the necessary flexibility for efficient implementation of concatena- 
tion and splitting, as we shall see later. 

Adding a new element to the left of the sequence is trivial, except when the initial 
buffer already contains four elements. In that case, we push three of the elements 
as a Node, leaving two behind: 



The function '<' requires a non-schematic form of recursion, called polymorphic 



3.2 Deque operations 



infixr 5 < 




a — ► FingerTree a — > FingerTree a 
Single a 

Deep [a] Empty [b] 

Deep [a, b] (Node3 c d e <m) sf 

Deep ([a] -H- pr) m sf 
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recursion (Mycroft, 1984): the recursive call adds a node (not an element) to a 
finger tree that contains Nodes. 

Adding to the right end of the sequence is the mirror image of the above: 

infixl 5 > 

(>) :: FingerTree a — > a — > FingerTree a 

Empty > a = Single a 

Single b > a — Deep [b] Empty [a] 

Deep pr m [e, d, c, b]> a = Deep pr (m > Node3 e d c) [b, a] 

Deep pr m sf > a = Deep pr m (sf 4f [a]) 

We shall also require the liftings of '<' and '>': 

(<') :: (Reduce f) =>■ / a — > FingerTree a — > FingerTree a 
(<') = reducer (<) 

(>') :: (Reduce f) =>■ FingerTree a — > / a — > FingerTree a 
(>') = reducel (>) 

Conversion to a finger tree may be defined using either '<" or '>"; we use the former: 

toTree :: (Reduce /) =>■ / a — » FingerTree a 
toTree s — s <f Empty 

To deconstruct a sequence, we define a type that expresses a view of the left end 
of a sequence as either empty or an element together with the rest of the sequence: 

data V^iewL s a = MZz, | Cons^ a (s a) 

The left view of a sequence maps a sequence to this type: 

viewL FingerTree a — > Viewi, FingerTree a 

viewL Empty = NUl 

vieuiL (Single x) = Cousl x Empty 

viewL (Deep pr m sf ) = Cons^ (head pr) (deep L (tail pr) m sf) 

Since the prefix pr of a Deep tree contains at least one element, we can get its head 
in constant time. However, the tail of the prefix may be empty, and thus unsuitable 
as a first argument to the Deep constructor. Hence we define a smart constructor 
that differs from Deep by allowing the prefix to contain zero to four elements, and 
in the empty case uses a left view of the middle tree to construct a tree of the 
correct shape: 

deep L :: [a] — ► FingerTree (Node a) — > Digit a — > FingerTree a 

deep L [] m sf = case view l m of 

Nil L toTree sf 

ConsL a m' — > Deep (toList a) m' sf 

deep L pr m sf — Deep pr m sf 

Then we can define separate functions isEmpty, headh, and tailh using the view: 
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isEmpty :: FingerTree a — ► Bool 

isEmpty x = case view l x of Nil l — > TWie 

Cons l — ► Fafee 

headr, :: FingerTree a — > a 

/lead l a; = case view i of Cons £, a _ — > a 

tai/i :: FingerTree a — > FingerTree a 

tail i, x = case view^ x of ConsL - x' — > a;' 



In a lazy language like Haskell, the tail part of the view is not computed unless it 
is used. In a strict language, it might be more useful to provide the three separate 
functions as primitives. 

We omit the code for the right view, which is a mirror image of the left view, 
including a function deep R corresponding to deep L . 

Each deque operation may recurse down the spine of the finger tree, and thus take 
O(logn) time in the worst case. However, it can be shown that these operations take 
only 6(1) amortized time, even in a persistent setting. The analysis is essentially 
the same as Okasaki's for implicit deques (Okasaki, 1998). We outline the main 
points here. 

We classify digits of two or three elements (which are isomorphic to elements 
of type Node a) as safe, and those of one or four elements as dangerous. A deque 
operation may only propagate to the next level from a dangerous digit, but in doing 
so it makes that digit safe, so that the next operation to reach that digit will not 
propagate. Thus, at most half of the operations descend one level, at most a quarter 
two levels, and so on. Consequently, in a sequence of operations, the average cost 
is constant. 

The same bounds hold in a persistent setting if subtrees are suspended using lazy 
evaluation. This ensures that transformations deep in the spine do not take place 
until a subsequent operation needs to go down that far. Because of the above prop- 
erties of safe and dangerous digits, by that time enough cheap shallow operations 
will have been performed to pay for this expensive evaluation. 

This argument may be formalized using Okasaki's debit analysis (1998) by as- 
signing the suspension of the middle tree in each Deep node as many debits as the 
node has safe digits. (It may have 0, 1 or 2.) A deque operation that descends k 
levels turns k dangerous digits into safe ones, and thus creates k debits to pay for 
the work involved. Applying Okasaki's technique of debit passing to any debits al- 
ready attached to these k nodes, one can show that each operation must discharge 
at most one debit, so the deque operations run in 0(1) amortized time. 

Lazy evaluation is the norm in a language such as Haskell. In a strict language 
that provides a lazy evaluation primitive, we need only suspend the middle subtree 
of each Deep node, so only ©(log n) suspensions are required in a tree of size n. Even 
in a lazy language, some space could be saved in practice by ensuring that these 
were the only suspensions in the tree, for example by using Haskell's strictness 
annotations. Indeed we might wish to make the operations even more strict, for 
example by forcing the middle subtree m in the recursive case of '<' to avoid building 
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a chain of suspensions. According to the above debit analysis, this suspension has 
been paid for by this time, so the amortized bounds would be unaffected. 

While we chose to push a NodeS in '<' and '>', leaving two elements behind, 
pushing a Node2 instead also yields the same bounds. Our choice means that the 
tree is mostly populated with nodes of degree 3, and therefore slightly shallower. 
It is also faster in queue-like applications, where elements are mostly added at one 
end and mostly removed from the other, because it means that only every third 
operation propagates to the next level. 

Exercise 1. In the above presentation, we represented digits as lists for simplicity. 
A more accurate and efficient implementation would use the type 

data Digit a = One a 
Two a a 
| Three a a a 
Four a a a a 

Rework the above definitions to use this definition of Digit. □ 

3.3 Concatenation 

We now consider concatenation of sequences represented as finger trees. The tree 
Empty will be an identity for concatenation, and concatenation with Singles reduces 
to '<' or '>', so the only difficult case is concatenation of two Deep trees. We can 
use the prefix of the first tree as the prefix of the result, and the suffix of the second 
tree as the suffix of the result. To combine the rest to make the new middle subtree, 
we require a function of type 

FingerTree {Node a) — ■> Digit a — ► Digit a — ► FingerTree (Node a) — ► 
Finger Tree (Node a) 

Since it is easy to combine the two digit arguments into a list of Nodes, we can 
obtain a recursive function by generalizing the concatenation function to take an 
additional list of elements: 

appS :: FingerTree a —* [a] — > FingerTree a — ► FingerTree a 

app3 Empty ts xs = ts <' xs 

app3 xs ts Empty = xs >' ts 

app3 (Single x) ts xs = x <(ts <' xs) 

app3 xs ts (Single x) = (xs >' ts) > x 

app3 (Deep pr l mi s/j) ts (Deep pr 2 m 2 sf 2 ) 

= Deep pr 1 (app3 m\ (nodes (sf \ -H- ts 4f pr 2 )) rn 2 ) sf 2 

where '<" and '>" are the previously defined liftings of '<' and '>', and nodes groups 
a list of elements into Nodes: 

nodes :: [a] — > [Node a] 

nodes [a, b] = \Node2 a b] 

nodes [a, b, c] = [NodeS a b c] 

nodes [a, b, c, d] = \Node2 a b, Node2 c d] 
nodes (a : b : c : xs) = NodeS a b c : nodes xs 
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Now we can concatenate two finger trees by starting with an empty list between 
them: 

(tx) :: FingerTree a — > FingerTree a — ► FingerTree a 
xs txl ys — appS xs [] ys 

It is easy to check that in each of the invocations of app5 arising from W, the list 
argument has length at most 4, so each of these invocations takes 0(1) time. The 
recursion terminates when we reach the bottom of the shallower tree, with up to 4 
insertions, so the total time taken is 0(log(min{ni, n 2 })), where n\ and n-i are the 
numbers of elements in the argument trees. 

Note that the debit analysis of Section 3.2 remains valid, because although W 
must discharge any debits on the nodes on the spines that are merged, there will 
be at most 0(log(min{ni, TI2})) of these. 

Exercise 2. The definition of nodes prefers to use NodeS where possible, but any 
other choices would still yield the same asymptotic bounds. If nodes were redefined 
to prefer Node2 , what would be the bound on the length of the list arguments in 
the invocations of app3 arising from W ? □ 
Exercise 3. Consider replacing the Digit type as discussed in Exercise 1, and also 
using a variant of app3 for each possible size of the middle argument. This leads to 
a faster but much larger definition. Hint: it is both simpler and less error-prone to 
write a program that generates the Haskell source for this definition. □ 

4 Annotated sequences 

Finger trees serve remarkably well as implementations of catenable deques. In this 
section we shall present a modification of these trees, which also provides positional 
and set-theoretic operations, such as taking or dropping the first n elements of a 
sequence or partitioning an ordered sequence around a pivot element. 

As a common characteristic, all of these operations involve searching for an el- 
ement with a certain property. For example, partitioning an ordered sequence in- 
volves searching for the minimum element that is larger than or equal to the pivot 
element. Thus, if we want to implement the additional operations with reasonable 
efficiency, we must be able to steer this search. To this end we augment the internal 
nodes of a finger tree by an additional field that contains positional or ordering 
information or both. 

4-1 Measurements 

We take the fairly liberal view that an annotation represents a cached reduction 
with some monoid. Reductions, possibly cached, are captured by the following class 
declaration (for clarity, we use mixfix notation, which is not possible in Haskell). 

class (Monoid v) => Measured a v where 
|H|::a-« ' 

Think of a as the type of some tree and v as the type of an associated measure- 
ment. As customary, v must possess a monoidal structure so that we can easily 
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combine measurements of subtrees independent of their nesting. For instance, the 
'size' measure maps onto the monoid of natural numbers with addition. In this par- 
ticular example, the binary operation is not only associative but also commutative, 
but this need not be the case in general. 

The class Measured is parameterized by both a and v and represents a true many- 
to-many relation on types: the same measurement may be imposed on a variety of 
types; the same type may be measured in many different ways. Unfortunately, this 
means that the class declaration is not legal Haskell 98, which only admits single- 
parameter type classes. Multiple parameter classes (Peyton Jones et at, 1997), 
however, are supported as an extension by the most widely used Haskell implemen- 
tations. 



4-2 Caching measurements 

Obtaining a measure should be cheap. We shall ensure that finger trees support 
this operation with a bounded number of '©' operations. To this end we cache 
measurements in 2-3 nodes. 

data Node v a = Node2 v a a \ NodeS v a a a 

The smart constructors node2 and nodeS automatically add an explicit annotation 

node2 :: (Measured a v) => a — > a — > Node v a 

node2 a b = Node2 (\\a\\ © \\b\\) a b 

node3 :: (Measured av)^>a^>a^a^> Node v a 

node3 a b c = NodeS (\\a\\ © ||&|| © ||c||) a b c 

This annotation is then retrieved by the measure function: 

instance (Monoid v) Measured (Node v a) v where 
\\Node2 v _ _|| = v 
\\Node3 v || = v 

The smart constructors and the instance declaration are completely generic: they 
work for arbitrary annotations. 
Digits are measured on the fly. 

instance (Measured a v) =>■ Measured (Digit a) v where 
||.t,s|| = reducel (Xi a — > i © ||a||) 0 xs 

Recall that a digit is a buffer of one to four elements. Because the length of the 
buffer is bounded by a constant, the number of '©' operations is also bounded. 
Another possibility is to cache the measure of a digit, adding to the cost of digit 
construction but yielding a saving when computing the measure. The choice between 
these strategies would depend on the expected balance of query and modification 
operations, but they would differ only by a constant factor. 

Finger trees are modified in a similar manner to 2-3 nodes. Again there is a choice 
of caching strategies, but it suffices to cache deep nodes: 
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data FingerTree v a = Empty 
| Single a 

| Deep v (Digit a) (FingerTree v (Node v a)) (Digit a) 

The top level of an annotated finger tree contains elements of type a, the second 
level of type Node v a, the third of type Node v (Node v a), and so on. Whereas 
the tree type a changes from level to level, the annotation type v remains the same. 
In a nutshell, FingerTree is nested in a, but regular in v. 

Again, a smart constructor takes care of filling in the missing annotation 

deep :: (Measured a v) => 

Digit a — > FingerTree v (Node v a) — > Digit a — > FingerTree v a 
deep pr m sf = Deep (\\pr\\ © ||m|| © ||s/||) pr m sf 

This annotation is then available for subsequent measurements: 

instance (Measured a v) Measured (FingerTree v a) v where 
\\Empty\\ =0 
\\Singlex\\ =\\x\\ 
\\Deep v || = v 



4-3 Construction, deconstruction and concatenation 

None of the sequence operations introduced in Section 3 benefits from the anno- 
tations. Consequently, their re-implcmcntations can simply ignore the additional 
field in pattern matching and recompute it when constructing nodes. The latter 
is most easily accomplished by replacing the original constructors by their smart 
counterparts (in expressions). The re-implementation of '<' illustrates the necessary 
amendments: 
infixr 5 < 

(<) :: (Measured a v) => a — > FingerTree v a — > FingerTree v a 

a < Empty = Single a 

a< Single b = deep [a] Empty [b] 

a < Deep _ [6, c, d, e] m sf = deep [a, 6] (nodeS c d e < m) sf 

a < Deep _ pr m sf = deep ([a] 4f pr) m sf 

The remaining operations are changed accordingly. Their time bounds remain as 
in Section 3, except that they now refer to the number of '©' operations required 
to compute the new annotations. 



4-4 Splitting 

A scan is an iterated reduction that maps the sequence x\ , . . . , x n of elements to the 
sequence v\, . . . , v n of accumulated measurements where Vj = i © || © ■ • • © \\xj \\ 
for some start value i. The next operation we introduce will use the annotations 
to efficiently split a sequence based on properties of the scan with the monoid. 
For instance, a 'size' scan yields the positions of elements so that we can split a 
sequence at a given position. Since '©' is associative, a scan can be calculated on 
the fly using the cached reductions in the nodes. 
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To represent a container split around a distinguished element, with containers of 
elements to its left and right, we introduce the data type 

data Split f a = Split (/ a) a (/ a) 

The Split type singles out an element and its context, in a similar way to a finger 
in an imperative setting or a zipper in a functional setting. Splitting a finger tree 
provides us with a third finger besides the ones at the left and right ends of the 
tree. 

The function splitTree, which splits a finger tree, expects three arguments: a 
predicate p on measurements, an accumulator i, and a finger tree t. It is specified 

->p i A p(i® ||i||) => 

let Split I x r = splitTree p i t in toList I -ff [x] -ff toList r = tohist t 

A ->p(i® \\l\\) A p (i © ||/|| © ||.t||) 

The first conjunct states that Split I x r is actually a split of t, that is, the order 
and multiplicity of elements is preserved. As usual for this type of condition, it will 
be fairly obvious from the program text that the property is satisfied. Indeed, we 
take great care that the relative order of elements is preserved in the syntax of our 
code. 

The other conjunct states that splitTree splits t at a point where p applied 
to the accumulated measurement changes from False to True. The precondition 
expresses that this measurement is False at the left end and True at the right 
end, so there is at least one such point. The precondition implies, in particular, 
that t is non-empty Splitting is most useful if there is exactly one such point, 
but we need not assume this. Uniqueness of the split is guaranteed for monotonia 
predicates satisfying p x => p (x(By) or cquivalently ~^p (x © y) =>• ->p x for all x 
and y. Alternatively, an application might maintain an invariant on sequences that 
guarantees a unique split for predicates of interest; Section 4.7 provides an example 
of this. 

The picture below illustrates the functioning of splitTree. 
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where Vj = i © ||a;i|| © • • • © \\xj\\. In the diagram, Xi is the distinguished element 
where p changes from False to True. 

The helper function splitDigit, which satisfies an analogous specification, nicely 
illustrates splitting along a scan. 
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splitDigit :: (Measured a v) => (v — > Bool) — > u — > Digi£ a — > SpZzi [] a 
splitDigit p i [a] = Split [] a [] 
splitDigit p i (a : as) 

| p i' = Split [] a as 

| otherwise = let Split I x r = splitDigit p i' as in Split (a : I) x r 
where i' = i © II all 



The type signature makes precise that splitting a digit may produce empty left 
and right buffers ('[]' is the list type constructor). The reader is invited to check 
that the invariant expressed in the specification is indeed preserved through the 
recursive calls. Since splitDigit scans the list from left to right, it actually splits 
the list at the first point where p flips. The function splitTree proceeds in a similar 
fashion, but it may skip whole subtrees (which become deeper as we descend the 
tree). This behaviour is necessary to make splitTree run in logarithmic time, but 
it also means that there is no guarantee that the calculated split is the first one. 
The definition of splitTree uses the auxiliary function deep L from Section 3.2 and 
its mirror image deep R . 



splitTree :: (Measured a v) => 

(v — > Bool) — > v — > FingerTree v a — > Split (FingerTree v) a 
splitTree p i (Single x) — Split Empty x Empty 
splitTree p i (Deep _ pr m sf) 

| p vpr = let Split I x r = splitDigit p i pr 

in Split (toTree I) x (deep L r m sf) 
| p vm = let Split ml xs mr = splitTree p vpr m 

Split I x r = splitDigit p (vpr © \\ml\\) (toList xs) 
in Split (deep R pr ml I) x (deep L r mr sf) 
| otherwise = let Split I x r = splitDigit p vm sf 

in Split (deep R pr ml) x (toTree r) 
where vpr = i © ||pr|| 
vm — vpr © || m || 



The cost of splitTree is proportional to the depth of the split point in the input 
finger tree, which as we observed in Section 3 is 8(logrf), where d is the distance to 
the nearest end. Thus the cost of splitTree is 0(log(min{n/, n r })) '©' operations, 
where n/ and n r are the sizes of the subtrees returned. As with W, this bound 
is also sufficient to discharge any debits on spine nodes that are split, so that the 
debit analysis of Section 3.2 remains valid. 

The above specification of splitTree is pleasingly simple, but we can strengthen it 
slightly. Plainly, splitTree requires a non-empty tree. An empty middle subtree m 
of a Deep node is not passed to a recursive call of splitTree, because in that case 
vm = vpr. Note also that the preconditions ->p i and p (i © ||f ||) arc used only if 
the split occurs at one or other end of the sequence. Consequently, if we treat those 
cases specially it suffices to require that the initial tree t is non-empty: 
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splitTree p i t in toList I -H- [a;] -ff toList r = toList t 
A (I = Empty V ->p (i © ||f||)) 
A (r = Empty V p (i ® \\l\\ © ||a;||)) 

Thus the trees to the left and to the right are either empty, in which case we can't 
say anything about p, or they satisfy the previous conditions. This specification 
allows us to define a function split that produces two sequences, specified by 

let (I, r) = split p xs in toList I 4f toList r = toList xs 

A (I = Empty V ->p \\l\\) 

A (r = Empty V p (\\l\\ © \\head L r||)) 

and implemented as follows: 

split :: (Measured a v) => 

(v —> Bool) —> FingerTree v a — > (FingerTree v a, FingerTree v a) 
split p Empty = (Empty, Empty) 
split p xs 

| p \\xs\\ = (I, x < r) 
| otherwise — (xs, Empty) 
where Split I x r = splitTree p 0 xs 

Note that we don't need any preconditions on p. If we relied on the weaker specifi- 
cation, we would require ~^p 0 even though p is never applied to 0 in the code. As 
an aside, the only monotonic predicate that holds for 0 is the one that is always 
true: when p is monotonic, p 0 implies p x for all x. 

Finally, we define convenient shortcuts for the left and right part of a split. 

take Until, drop Until :: (Measured a v) => 

(v — ► Bool) — > FingerTree v a — > FingerTree v a 
takeUntil p = fst ■ split p 
dropUntil p = snd ■ split p 

Exercise 4. Give a more efficient version of the following function: 

lookupTree :: (Measured a v) => (v — > Bool) — > v — > FingerTree v a — > (v, a) 
lookupTree p i t = let Split I x r = splitTree p i tin (i © ||/||, x) 

assuming ->p i A p (i © ||£||). □ 



t =/= Empty =>■ 

let SpZif / ir r 



^.5 Application: random- access sequences 

Let's move on to our first application. An efficient implementation of sequences 
ought to support fast positional operations such as accessing the nth element or 
splitting a sequence at a certain position. To this end we annotate finger trees with 
sizes. 
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newtype Size = Size{getSize 

deriving (Eq, Ord) 

instance Monoid Size where 
0 = Size 0 



Size m © Size n = Size (m + n) 

As announced in Section 2.1 we introduce a new type because the underlying type N 
of natural numbers is the carrier of several different monoids. A second new type 
is needed for the elements of the sequence. 

newtype Elem a = Elem{getElem :: a} 

Think of Elem as a marker for elements. Now, we can define 

newtype Seq a = Seq (FingerTree Size (Elem a)) 

instance Measured (Elem a) Size where 
\\Elem _|| = Size 1 

The instance declaration serves as a base case for measuring sizes: the size of an 
element is one. 

In a library, the types Size and Elem would be hidden from the user by defining 
suitable wrapper functions. Note in this respect that the Haskell Report guarantees 
that the use of newtypes does not incur any run-time penalty. 

The length of a sequence can now be computed in constant time. 

length :: Seq a — > N 

length (Seq xs) = getSize \\xs\\ 

Given these prerequisites, splitting a sequence at a certain position is simply a 
matter of calling split with the appropriate arguments. 

split At :: N — > Seq a — ► (Seq a, Seq a) 

splitAt i (Seq xs) = (Seq I, Seq r) 
where (I, r) = split (Size i <) xs 

As an aside, the predicate (Size i <) is monotonic as i < m implies i < m + n for 
natural numbers m and n, so in this case the split is unique. 
The implementation of indexing is equally straightforward. 



where Split _ x _ = splitTree (Size i <) (Size 0) xs 

Exercise 5. Use the definition of the Seq type to define a specialization of the deque 
operations, noting that it is sometimes faster to recompute sizes using subtraction 



A second interesting choice is to annotate with the maximum operation. Though 
this operation is associative, it might not have an identity. We can turn it into a 
monoid by adjoining an identity element, here MInfty (minus infinity). 




:: Seq a — ► N — > a 
= getElem x 



instead of addition. 



□ 



4-6 Application: max-priority queues 
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data Prio a = MInfty \ Prio a 

deriving (Eq, Ord) 

instance (Ord a) => Monoid (Prio a) where 
0 = MInfty 

MInfty ®p = p 

p © MInfty = p 

Prio m (B Prio n = Prio (m 'max 1 n) 

Unsurprisingly, measuring with this monoid gives us max-priority queues. 

newtype PQueue a = PQueue (FingerTree (Prio a) (Elem a)) 

instance (Ord a) => Measured (Elem a) (Prio a) where 
|| Elem x\\ — Prio x 

Interestingly, the type PQueue is quite close to 2-3-4 heaps (Cormen et at, 2001), 
which are priority queues based on 2-3-4 trees: all the data is stored in the leaves; 
the inner nodes hold the largest keys of their subtrees. As is typical of heaps, there 
is no constraint on the order of the elements within the leaves. Consequently, to 
insert an element we simply add it to the queue using '<' or '>'. Extracting the 
element with the largest key is implemented in terms of splitTree (we assume that 
the queue is non-empty). 

extractMax :: (Ord a) =^> PQueue a — > (a, PQueue a) 

extractMax (PQueue q) — (x, PQueue (I tx r)) 
where Split I (Elem x) r = splitTree (\\q\\ 0 q 

Exercise 6. Show that it is possible to find an element ^ any k in 9(log(n)) time, 
and all m of them in 6(mlog(n/m)) time. □ 

4-7 Application: ordered sequences 

The binary operation that always selects its second argument is also associative. As 
with max in the last section, it can be made into a monoid by adjoining an identity 
element, here called NoKey. 

data Key a = NoKey | Key a 

deriving (Eq, Ord) 

instance Monoid (Key a) where 
0 = NoKey 

k © NoKey = k 
_ffi k = k 

In a tree annotated using this monoid, each subtree is annotated with the last key 
it contains. If wc maintain the sequence in key order, we have an implementation 
of ordered sequences, with these annotations serving as split or signpost keys: 

newtype OrdSeq a = OrdSeq (FingerTree (Key a) (Elem a)) 

instance Measured (Elem a) (Key a) where 
\\Elem x\\ — Key x 
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The Key monoid can be seen as an optimization of the Prio monoid: if the sequence 
is ordered, maximum simplifies to 'second' with NoKey playing the role of MInfty. 

Ordered sequences subsume priority queues, as we have immediate access to the 
smallest and the greatest element, and search trees, as we can partition ordered 
sequences in logarithmic time. 

The first operation on ordered sequences partitions before the first element ^ k: 

partition :: (Ord a) => a — > OrdSeq a — > (OrdSeq a, OrdSeq a) 
partition k (OrdSeq xs) = (OrdSeq I, OrdSeq r) 
where (I, r) = split Key k) xs 

Perhaps surprisingly, insertion and deletion can also be defined in terms of split. 
The implementations below are both logarithmic, though special purpose versions 
would have smaller constant factors. 

insert :: (Ord a) => a — > OrdSeq a — > OrdSeq a 
insert x (OrdSeq xs) = OrdSeq (I m (Elem x<r)) 
where (I, r) = split Key x) xs 

deleteAll :: (Ord a) => a — > OrdSeq a — > OrdSeq a 
deleteAll x (OrdSeq xs) = OrdSeq (I ix r') 
where (I, r) — split Key x) xs 
(_, r') = split (> Key x) r 

The non-strict split of xs yields a third finger pointing to the position where x must 
be inserted or deleted. This technique is reminiscent of splay trees (Sleator & Tarjan, 
1985) except that the subsequent append (W) does away with the additional finger. 
Deletion is defined to remove an arbitrary number of elements. The cost of the first 
split in deleteAll is logarithmic in the distance to the nearest end, while the cost of 
the second is logarithmic in the number of matches. 

Building upon concatenation and split we can also define the merge of two ordered 
sequences. 

merge :: (Ord a) => OrdSeq a — ► OrdSeq a — ► OrdSeq a 
merge (OrdSeq xs) (OrdSeq ys) = OrdSeq (merge' xs ys) 
where merge 1 as bs = case view l bs of 

Nil L -» as 

ConsL a bs' — ► / x (a < merge' bs' r) 
where (/, r) = split (> ||a||) as 

The function partitions the two input sequences into the minimum number of seg- 
ments that must be reordered to obtain an ordered output sequence. The algorithm 
is similar to the merge of lists. It is, however, more efficient since it uses exponential 
and binary search rather than linear search. It can be shown that merge is asymp- 
totically optimal: it takes 6(mlog(n/m)) amortized time where m and n are the 
lengths of the shorter and longer input sequences respectively (Moffat et at, 1992). 
Furthermore, merge exhibits an adaptive behaviour: for sequences xs and ys with 
max xs ^ min ys (or vice versa) merging degenerates to concatenation. So, in this 
favourable case we have a running time of 9(logm). 
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Exercise 7. In the above implementation, all the annotations stored in the tree 
have the form Key x, which suggests an optimization. Redefine splitTree, deep and 
the Measured instance for lists to avoid relying on 0 being the right identity of '©', 
so we can use the binary function Xx y — > y instead. □ 

Exercise 8. Implement the intersection of two ordered sequences. □ 

Exercise 9. Combine sizes with split keys to get order statistics. The ith order 
statistic of a collection of n elements is the ith smallest element (Cormen et at, 
2001). □ 



4-8 Application: interval trees 

Our final application augments finger trees to support query operations on sets 
of intervals: we show how to find an interval that overlaps with a given interval in 
0(log(n)), and all m of them in 0(to log(n/m)) time. This application is particularly 
interesting as it makes use of two annotations that are maintained in parallel (a 
similar application is the subject of Exercise 9). For this we use the product of 
monoids, given by 

instance (Monoid a, Monoid b) Monoid (a, b) where 
0 =(0,0) 
(a, b) © (a', b') = (a © a', b © b') 

A closed interval i is an ordered pair of real numbers 

data Interval — Interval^ low :: R, high :: R} 

with low i < high i. Two intervals i and j overlap when 

low i ^ high j A low j < high i. 

An interval tree is a finger tree augmented by split keys and priorities. 

newtype IntervalTree = IntervalTree (FingerTree (Key R, Prio R) Interval) 

instance Measured Interval (Key R, Prio R) where 
= (Key (low i),Prio (high i)) 

We order the intervals by their low endpoints, and additionally annotate with the 
maximum of the high endpoints. Interval trees are an instance of priority search 
queues (McCreight, 1985), which simultaneously support dictionary operations on 
keys and priority queue operations on priorities. 

The predicates below are useful for querying keys and priorities. 

atleast, greater :: R — > (Key R, Prio R) — ► Bool 
atleast k (_, n) — Prio k ^ n 
greater k (n, _) = n > Key k 

The function intervals earch finds an interval overlapping with a given interval in 
the above times, if such an interval exists. 
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intervals earch :: IntervalTree —> Interval — > Maybe Interval 

intervals earch {IntervalTree t) i 

| atleast (low i) \\t\\ A low x ^ /wg/i i = Jwsi x 

| otherwise = Nothing 

where Split _ x _ = splitTree (atleast (low ij) 0 £ 

The splitTree gives us the interval x with the smallest low cndpoint whose high 
endpoint is at least the low endpoint of the query interval — the initial conjunct 
atleast (low i) \\t\\ ensures that such an interval exists. It then remains to check that 
low x ^ high i. If this test fails, then we know that there can't be any intersecting 
interval: any x' to the left of x in the sequence doesn't satisfy low i < high x'; any 
x' to the right doesn't satisfy low x' ^ high i. 
Finding all intervals is just as easy: 

intervalMatch :: IntervalTree — ► Interval — > [Interval] 

intervalMatch (IntervalTree t) i — matches (takeUntil (greater (high ij) t) 
where matches xs = case viewL (drop Until (atleast (low i j) xs) of 
Nil L -» [] 

ConsL x xs' — > x : matches xs' 

The application of takeUntil selects the intervals x with low x ^ /wg/i i, so it 
remains to select the intervals x from this sequence that also satisfy high x < low i. 
However, these intervals are not necessarily adjacent in the sequence, because it 
is ordered by the low endpoint of the intervals. For example, a long interval that 
extends to i might be followed in the sequence by a shorter one that does not. 
Hence we use dropUntil to find the next intersecting interval. 



5 Related work and conclusion 

Finger trees were originally defined by reversing certain pointers in a search tree to 
accelerate operations in the vicinity of distinguished positions in the tree, known 
as fingers (Guibas et at, 1977). The original versions used various forms of B- 
trees, including 2-3 trees, variants using AVL trees have subsequently appeared 
(Tsakalidis, 1985). 

Our functional 2-3 finger trees are an instance of a general design technique in- 
troduced by Okasaki (1998), called implicit recursive slowdown. We have already 
noted that these trees are an extension of his implicit deque structure, replacing 
pairs with 2-3 nodes to provide the flexibility required for efficient concatenation 
and splitting. Okasaki also presents a wealth of functional implementations of se- 
quences and other data types. Implementations of functional deques with constant 
time concatenation are known (Kaplan & Tarjan, 1995; Okasaki, 1997), but these 
structures do not admit efficient insertion, deletion or splitting. Kaplan and Tarjan 
(1996) defined two functional data structures with the same time bounds as here, 
except that their bounds are worst case, and do not require laziness. They achieved 
this using structures based on segmented number systems, which are considerably 
more involved than the structure used here. Their second implementation is iso- 
morphic to a restricted version of 2-3 finger trees, though their representation using 
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Tabic 1. Comparing persistent sequence implementations 



segmented spines is vastly more complex, as are the operations required to preserve 
their invariants. They also sketched a third implementation, claimed to provide a 
doubly logarithmic concatenation time, but many details of this variant are unclear. 

The biggest advantage of 2-3 finger trees is the simplicity of the data structure 
and operations: the implementations of the deque operations follow naturally from 
the structure of finger trees, and so do the implementations of concatenation and 
splitting. In addition, first experiences with a prototypical implementation suggest 
that finger trees perform well in practice. 

For example, Table 1 compares times for a number of persistent implementa- 
tions of sequences. The times are averaged over 5,000,000 operations on a sequence 
of 10,000 elements, on a Pentium III 450 with GHC 6.4. Concatenation is not 
considered, as it is difficult to separate the cost of building a sequence using con- 
catenations from that of accessing it. We conclude that finger trees offer a general 
purpose implementation of sequences that is competitive with existing persistent 
implementations of subsets of the operations, though ephemeral implementations 
will usually be faster. As search trees, finger trees seem to be 3 to 5 times slower 
than the best balanced binary tree implementations. Nevertheless, they may be 
competitive when a wider set of operations is required. 

The treatment of annotations nicely demonstrates the benefits of abstraction: a 
single generic split function serves a multitude of purposes. 
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